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ABSTRACT 



In this paper we derive the galaxy luminosity function from the Nearby Optical Galaxy (NOG) 
sample, which is a nearly complete, magnitude-limited (B <14 mag), all-sky sample of nearby 
optical galaxies (~6400 galaxies with cz < 5500 km s _1 ). For this local sample, we use galaxy 
distance estimates based on different peculiar velocity models. Therefore, the derivation of 
the luminosity function is carried out using the locations of field and grouped galaxies in real 
distance space. 

The local field galaxy luminosity function in the B system is well described by a Schechter 
function with a slope a ~ -1.1, a low normalization factor (<J>* ~ 0.006 Mpc~ 3 ), and a par- 
ticularly bright characteristic magnitude (Mjg ~ —20.6) (Hq = 75 kmT x Mpc~ v ). The exact 
values of the Schechter parameters slightly depend on the adopted peculiar velocity field models. 
Peculiar motion effects are of the order of statistical errors and cause at most variations of 0.8 
in a and 0.2 mag in M|j. Our MJ>-value is brighter by a few tenths of magnitude than previous 
corresponding values, because, referring to total corrected blue magnitudes, better represent the 
galaxy light. 

Also the selection function, evaluated in terms of the luminosity function, appears to be little 
sensitive to the adopted peculiar velocity field models, which, however, bear a large impact on 
the local galaxy density on the smallest scales. 

The shape of the luminosity function of spiral galaxies does not differ significantly from that of 
E-S0 galaxies. On the other hand, the late-type spirals and irregulars have a very steeply rising 
luminosity function towards the faint end (a ~-2.3 - -2.4), whereas the ellipticals appreciably 
decrease in number towards low luminosities. 

The presence of galaxy systems in the NOG sample does not affect significantly the field galaxy 
luminosity function, since environmental effects on the total luminosity function appear to be 
marginal. The luminosity function of the members of the richest galaxy systems tends to show 
a slightly brighter M^-value than the norm. 

In the light of constraints imposed by the observed galaxy number counts, the low normaliza- 
tion of the luminosity function suggests that the nearby universe (cz <j5000 km/s) examined in 
this paper may be underdense by a factor ~1.5. 

Subject headings: galaxies: distances and redshifts - galaxies: fundamental parameters - 
galaxies: luminosity function, mass function - cosmology: observations 
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1 INTRODUCTION 



Much work has gone into the determination of the local (z = 0) optical galaxy luminosity 
function (LF) (see, e.g., the review by Binggeli, Sandage & Tammann 1988 for early works), 
which is a fundamental quantity in cosmology. However, there is still uncertainty in its normal- 
ization and in its faint end. These quantities are important for the question of whether there is 
significant evolution in galaxy properties to z ~ 0.5 and for the explanation of the faint galaxy 
counts. As stressed in the reviews by Colless (1997) and Ellis (1997), many efforts are still 
needed to verify the detailed form of the local LF. Furthermore, in a magnitude-limited galaxy 
sample the knowledge of the LF is an essential ingredient for evaluating the selection function 
for the sample. 

In this paper we determine the optical galaxy LF — and, hence the intimately related selection 
function — of a magnitude-limited complete sample of nearby optical galaxies. The present work 
is indeed meant to be an important step towards the recovering of the optical galaxy density field 
in the nearby universe (see, e.g., Marinoni et al. 1998a for a preliminary account of this work). 
As previously discussed (e.g., Hudson 1993, Santiago et al. 1995), optical galaxy samples are 
more suitable for mapping the galaxy density field on small scales than IRAS-selected samples, 
which have been frequently used as tracers of the galaxy density field on large scales. 

In this work we consider the magnitude-limited, all-sky sample of nearby galaxies which was 
extracted by Garcia et al. (1993) from the Lyon-Meudon Extragalactic Database (LEDA). 
This sample comprises 6392 galaxies with recession velocities cz <5500 km s _1 and corrected 
total blue magnitudes B <14 mag. Although the different optical galaxy catalogues, from 
which data are collected and homogenized in the LEDA, have different limits of completeness 
in apparent magnitude or angular diameter, the above-mentioned galaxy sample was found to 
be substantially complete down to its limiting magnitude B = 14 mag, for galactic latitudes 
|6| > 20°) (Garcia et al. 1993). On the basis of the LEDA compilation of bright galaxies with 
and without redshifts, we estimate that the completeness level of this sample limited to |6| > 20° 
(5366 galaxies) is ~85% . For the sample of 5832 galaxies having \b\ > 15°, the completeness 
level decreases to ~80%, according to the deficit of galaxy counts at low |6|-values. Most the 
bright galaxies with unknown redshift lie at |6| < 30°. 

Group assignments for the galaxies of this sample have been already provided by Garcia (1993), 
who identified groups by employing the percolation or friends of friends method proposed by 
Huchra & Geller (1982) and the hierarchical clustering method (e.g., Tully 1987). The adopted 
final catalog of groups was defined as that one which includes only groups common to the two 
catalogs. Accordingly, we shall treat 3381 ungrouped galaxies as field galaxies. Of the 3011 
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members of the 485 groups (with at least three members), in the present study we omit 22 
galaxies which are possible members (members which almost satisfy the selection criteria). 

Garcia et al. (1993) tabulated several parameters for each galaxy, such as the morphological 
type, the corrected total blue magnitude B and the corrected angular size (the latters trans- 
formed in the standard systems of the RC3 catalog by de Vaucouleurs et al. 1991), and the 
distance modulus. In general, this quantity has been derived simply from redshift-dependent 
distances with an adopted Hubble constant Hq = 75 km s^^-Mpc" 1 , except for the nearest 
galaxies (with cz <1500 km s _1 ), where redshift-independent distances obtained from some 
distance indicators (DIs) (i.e. from a blue Tully-Fisher relation which relates luminosities to 
maximum rotation velocities and a luminosity-luminosity index relation) have been also par- 
tially considered. 

In Paper I (Marinoni et al. 1998b), we corrected the redshift-distances of the galaxies and 
groups of this sample for peculiar motions by means of models of the peculiar velocity field. 
We purposely avoided correcting distances by calculating the peculiar velocity field on the basis 
of the positions and redshifts of the galaxy sample itself (in the linear theory approximation), 
because the galaxy sample does not presumably include all relevant gravitational sources (such 
as the Shapley concentration) for local peculiar motions. We instead followed two basic inde- 
pendent approaches to modeling the velocity field in the nearby universe: i) a modified cluster 
dipole model, which is Branchini & Plionis' (1996) optical cluster dipole reconstruction scheme 
modified by the inclusion of a local model of the Virgocentric infall in the Local Supercluster 
region, ii) a semi-linear approach which uses a multi-attractor model fitted to the Mark II and 
Mark III catalogs of galaxy peculiar velocity (Willick et al. 1997). In this model the veloc- 
ity field was assumed to be generated by a few prominent gravitational sources (Virgo cluster, 
Great Attractor, Perseus-Pisces and Shapley superclusters) , which were described as spherically 
symmetric masses with King-type mass density radial distributions. 

Furthermore, we inverted the redshift-distance relations relative to different velocity field 
models and solved the problem of the triple- valued zones of this relation (these zones can appear 
in the vicinity of prominent overdensities) by using blue Tully-Fisher relations calibrated on 
suitably defined samples of objects having distances predicted by peculiar velocity models. 

In this way, in Paper I we provided homogeneous estimates of distances for the individual 3381 
field galaxies and the 485 groups of the above-mentioned galaxy sample, hereinafter denoted as 
Nearby Optical Galaxy (NOG) sample. 

There are three main conceptual differences between our present analysis of the local galaxy 
LF and previous relevant studies based on other, comparably shallow samples, such as the CfAl 
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(Davis & Huchra 1982), CfA2 (Marzke et al. 1994a), SSRS2 (da Costa et al. 1994), and ORS 
(Santiago et al. 1995, 1996). 

First, we use the distance estimates given in Paper I for determining the galaxy LF of the 
NOG sample. Thus, our derivation of the galaxy LF is carried out using the locations of field and 
grouped objects in real distance space, not in redshift space as is usually done in the literature. 
Having calculated different models of galaxy distances, we are also able to quantify to what 
extent differences in the current views on the cosmic flows can affect the determination of the 
galaxy LF of an optical galaxy sample. We will find that the luminosity and selection function 
of our optical and local sample are little sensitive to the adopted peculiar velocity field, thus 
confirming Yahil et al.'s (1991) finding based on IRAS galaxy samples (and different velocity 
field models). 

Second, we rely on apparent total B magnitudes fully corrected for internal extinction, Galactic 
extinction and K-dimming (Garcia et al. 1993). These corrections lead to bright magnitudes. In 
particular, the first correction, which is conspicuous in very inclined spiral galaxies, is generally 
neglected in generic redshift surveys, which comprise many faint galaxies of unknown inclination. 
The average correction for generic samples of bright spirals is ~0.2-0.3 mag in the B band (de 
Vaucouleurs, de Vaucouleurs Sz Corwin 1976; Hasegawa & Umemura 1993). Being selected on 
the basis of homogenized corrected magnitudes, NOG is in principle designed to provide a good 
estimate of the bright end of the local LF. 

Third, the NOG sample, which contains a large catalog of groups selected in a homogeneous 
way according to well-defined selection criteria, allows us to investigate on possible environmental 
effects on the LF, specifically on differences between the LFs of non-grouped (field) and grouped 
galaxies. 

A different, complementary approach to the definition of an all-sky sample of nearby optical 
galaxies with good completeness in redshift, the "Optical Redshift Survey" (ORS), was followed 
by Santiago et al. (1995). The ORS sample contains 8286 galaxies with known redshift and 
consists of two overlapping optically-selected samples (limited in apparent magnitude and diam- 
eter, respectively) which cover almost all the sky with |6| > 20°. Each sample is a concatenation 
of three subsamples drawn from the Uppsala General Catalogue (UGC) in the north, the Euro- 
pean Southern Observatory (ESO) catalogue in the south, and the Extension to the Southern 
Observatory Catalogue (ESGC) in a strip just south of the celestial equator. The authors se- 
lected their own galaxy sample according to the raw (observed) magnitudes and diameters and 
then quantified the effects of Galactic extinction (as well as random and systematic errors) on 
the resulting galaxy density field, which was calculated out to cz=8000 km/s in redshift space 
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(Santiago et al. 1996). Adding the IRAS 1.2 Jy galaxy sample (Fisher et al. 1995) in the 
unsurveyed zone of avoidance (|6| > 20°) and at large distances (cz >8000 km/s), Baker et al. 
(1998) calculated the resulting peculiar velocity field. 

Notwithstanding the differences in selection criteria and the larger redshift incompleteness of 
the NOG, the distribution of NOG galaxies on the sky is qualitatively similar to that of ORS 
galaxies; both samples delineate similar major structures in the nearby universe. 

The outline of our paper is as follows. In §2 we address the determination of the Schechter- 
type blue LF for the NOG sample. In §3 we describe our results concerning the total LF, the 
morphological type-specific LFs, and the comparison between the field and grouped galaxy LFs. 
In §4 we compare our results with previous LF determinations and discuss some implications 
related to galaxy number counts. In §5 we define the selection function of the NOG sample in 
terms of the galaxy LF. Conclusions are drawn in §6. 

Throughout, the Hubble constant is 75 /175 km s _1 Mpc^ 1 . 

2 The galaxy luminosity function: analysis 

We evaluate the galaxy LF of the NOG sample in the standard B system of the RC3 catalog. 
We rely on the various sets of distances (corrected for peculiar motions) given by Marinoni et al. 
(1998b) for field and grouped galaxies and on the total B magnitudes as derived by Garcia et al. 
(1993), who transformed the original raw data to the standard B system of the RC3 catalog. For 
numerous NOG galaxies total B magnitudes were carefully derived from aperture photometry 
and detailed (photographic or CCD) surface photometry. Garcia et al. (1993) corrected the 
total B magnitudes for internal extinction (following Tully & Fouque 1985), Galactic extinction 
(following Burstein & Heiles 1978 a,b, 1982, 1984), and K-dimming (Pence 1976). 

In general, extensive, bright magnitude-limited samples such as the NOG sample, which is 
currently one of the largest samples used in the LF determination, can provide good determi- 
nations of the local LF, except for its faint end (i.e., roughly for Mb > —15). Compared to 
other wide-angle, comparably shallow, bright magnitude-limited galaxy samples, such as the 
CfAl (Davis & Huchra 1982), CfA2 (Marzke et al. 1994a), and SSRS2 (da Costa et al. 1994; 
Marzke & da Costa 1997) redshift surveys, the NOG sample may be less sensitive to local density 
fluctuations because it covers a much larger solid angle. 

Deeper magnitude-limited redshift surveys are in principle better designed to address issues 
concerning the faint end of the local LF and its normalization, which in shallow galaxy samples 
could be affected by local density fluctuations. But, in practice, surveys which extend local 
samples to fainter magnitudes and higher redshifts may begin suffering more seriously from 



6 



some kinds of uncertainties: i) their photometric quality gets worse and, in particular, the 
photographic photometry scale may be not well-defined (see, e.g., Takamiya, Kron & Kron 
1995, Metcalfe, Fong & Shanks 1995, Rousseau et al. 1996, for systematic errors in Zwicky 
magnitudes, APM and COSMOS bj photometry, respectively); ii) their completeness gets more 
doubtful, because, at a given magnitude, they tend to miss more compact and low-surface 
brightness (LSB) galaxies by virtue of selection effects inherent in standard image detection 
algorithms (e.g., Disney & Philipps 1985; Davies, Disney & Phillipps 1989). As a matter of fact, 
at a given magnitude, galaxies which lie at nearer distances are selected over a broader interval 
of surface brightness (e.g., Phillipps, Davies, & Disney 1990). 

Much of the meaningful discrepancies between the local LFs derived in the literature (see 
§4 below) can be ascribed to photometric problems (magnitude scale errors, use of different 
isophotal or aperture magnitude systems, lack of adequate corrections of observed magnitudes) 
and to different observational selection criteria used in constructing galaxy samples. Both kinds 
of problems are a little minimized in NOG, because this sample comprises both bright and 
nearby galaxies. 

Furthermore, the deepest surveys suffer from redshift dependence in the measured isophotal 
or aperture magnitudes, due to the combined effects of the point spread function and surface 
brightness variation (e.g., Dalcanton 1998). In any case, very faint magnitude- limited redshift 
surveys (e.g., Ellis 1997), which are very important tools for probing evolutionary effects on the 
galaxy LF, have difficulties in giving unambiguous constraints on the local galaxy LF, because 
they are sensitive to evolutionary effects (e.g., the problem of the faint blue galaxy excess) and 
to the local normalization problem (e.g., Driver & Phillipps 1996). 

2.1 Estimating the shape 

We evaluate the shape of the galaxy LF following Turner's (1979) method (see also de Lap- 
parent, Geller & Hucra 1989). As recently demonstrated by Willmer (1997), who compared 
the robustness of different Schechter-type LF estimators by means of Monte Carlo simulations, 
this method — as well as all widely used methods, except the 1/V max technique — gives a 
substantially unbiased estimate of the shape of this function, especially in the case of large 
samples. 

Under the assumption that the LF is universal (i.e. independent of location in the universe), 
the galaxy number per unit volume and magnitude has a form separable into a product of 
functions of absolute magnitude M and position <I>(r, M)dMdr 3 = ip(M)p(r)dMdr 3 . 

In this notation, ip represents the fraction of objects with luminosity M, i.e. the shape of LF. 
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For a magnitude-limited sample of galaxies, the mean number of objects with luminosity M is 



n{M)dM 



p(r)dV 

V ma x(M) 



ip(M)dM = S(M)ip(M)dM (1) 



where V max (M) is the maximum volume in which galaxies brighter than M are still visible. 
Considering the number density of galaxies brighter than M 

rM 

N(< M) = S(M) / ip(M')dM' = S{M)V(M) (2) 

J — oo 

and taking the ratio between eq. [I] and eq. we are left with a quantity Y{M)dM which is 
density-independent. Rewriting this quantity using the following differential form 

Y(M)dM = -j^L (3) 
an analitic expression for the shape of the LF is given by 

rM 

cp(M) = (p Y(M)exp / Y(M')dM' (4) 
Jm 

where Mq is the brightest absolute magnitude we use in the estimation of the shape of the LF. 

Inside V max (M), and for small absolute magnitude increments (i.e. for small distance incre- 
ments 6r), Y(M) can be estimated using the following approximation (de Lapparent, Geller & 
Huchra, 1989) 

Y{M) W N{< M)5M (5) 

where iV(< M) is the number of galaxies brighter than M. Even if this estimation technique 
is intrinsically non parametric, in order to make comparisons with previous works, and also in 
order to avoid biases in the determination of the LF shape (de Lapparent et al., 1989), we fit 
eq. H against the one predicted, assuming, for the LF shape ip in the N(<M) estimator (eq. |2|), 
the Schechter (1976) form 

V? oc i 0- 4 (^ , '-^)(°+ 1 )exp[-10 - 4 ( M *- M )]. (6) 

In this equation, M* is a fiducial magnitude which characterizes the point (the "knee" of the 
function) where the form of the function changes from a power-law (with slope a) to an expo- 
nential slope 

With this least-square algorithm, weighted assuming Poisson fluctuations, we can measure the 
best parameters a and M* for the Schechter function. 
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However, in order to account for random errors in the magnitudes, we convolve the Schechter 
function $ s (ilf) with a Gaussian magnitude error distribution with zero mean and dispersion 
a m = 0.2 mag to give an observed LF <J> Q (M): 



1 r+oo 

* (M) = -== / * a (M') exp 



< M '- M l' d M> (7) 



2ai 



We perform the least-square fit of the LF estimator Y using various distance intervals of 
5r =250, 500 and 700 km/s. The lower limit for 5r was fixed by the spread caused by random 
noise in the peculiar velocity field, while the upper limit was set at the peak value of the absolute 
strength of the peculiar velocities. A variation in Sr over the range 250-700 km s _1 causes changes 
in the LF parameters a and M* which are of the order of uncertainties. 

2.2 Estimating the normalization 

Involving ratios between the differential and the integrated LFs, the estimator Y does not 
provide any information about the normalization factor 0* (i.e., the scaled number density of 
galaxies with luminosity M*) which has to be derived independently. We normalize the Schechter 
function using the fact that the mean spatial galaxy density n is related to the LF by the relation 

" ■ (8) 



Ms 

3>(M)dM 



where $ (M) is the convolved Schechter function with (f>* set equal to one. Since the LF is 
poorly constrained at the faint end, we cut off the integral at the limit of L s = L m i n (r s ), where 
Lmin(r) is the minimum luminosity necessary for a galaxy at distance r (in Mpc) to make it 
into the sample and r s is taken equal to r s = 500/(75 • /175) ~ 6-7hj§ Mpc. L m i n (r) corresponds 
to the absolute magnitude Mb = — 51ogr — 25 + Bn m , where Bn m = 14 mag is the limiting 
apparent magnitude of our galaxy sample; thus L s corresponds to M s = —15.12 + hloghj^. 

Various methods have been proposed to estimate the mean density. Given our complete and 
magnitude-limited sample of galaxy with r m j n < V{ < r max , we calculate the mean galaxy 
density n using the unbiased estimator 

Ngal 

f n = rr ma x 7 aiy \ ' (9) 

/ — I drS(r)w(r) 

Jr min \dr J 



where S(r) is the selection function given in equation |ll| below, f=0.8 is the adopted com 
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pleteness level (see §1) and w(r) is the weight function 



w(r) 



1 



(10) 



[l+4irnJ 3 (r c )F(r)Y 



,r 



which minimizes the variance in the estimate of n (Davis &: Huchra, 1982). For the second 
moment of the spatial two-point galaxy correlation function £(r), we have adopted the value 
J 3 = 320h~ 3 Mpc 3 obtained for r c = 40 h~ x Mpc (Marzke et al. 1994a). Even if the value of 
J3 is uncertain, because the contribution for r > r c could be substantial, the final result for n 
(and (j)*) depends little on the value of J3 used (Loveday et al. 1992). In our case, doubling or 
halving the adopted value of J3 makes a few per cent difference to n. 

In this minimum-variance weighting scheme, n is determined iteratively using eqs. || and 
|l0| while, for errors estimate, we add in quadrature the uncertainties in n arising from galaxy 
clustering (5n ~ n[4-7r J3/I/] ' 5 ) with those arising from varying the parameters of the luminosity 
function along their joint 1 a error ellipse (Lin et al. 1996). 

We have checked that the incompleteness of NOG does not appreciably depend on galaxy mor- 
phological type and mean surface brightness (inside the isophote 25 B-mag arcsec -2 ), whereas it 
gets weaker at brighter B magnitudes (e.g., /=0.90 (0.92) for |b| > 15° (|6| > 20°) and B <13.5; 
/=0.93 (0.96) for \b\ > 15° (\b\ > 20°) and B <13.0). After having derived in detail the degree 
of completeness of NOG as a function of the B magnitude, under the simple assumption of a 
uniform distribution of the missing galaxies in redshift space, we have calculated that this de- 
pendence does not appreciably affect the LF shape; it translates into a bias towards M*-values 
systematically brighter by 0.1 mag (and to a- values systematically flatter by 0.01 at most). 

In order to take into account border effects which arise from the mapping between redshift 
space and true-distance space and which influence the completeness of the sample (see Fig. 2 
below), we evaluate the Schechter-type LF for the galaxies with distances r < 5250/(75 /175) Mpc 
only. Further we consider only galaxies having \b\ > 15°, in order to avoid a large incompleteness, 
and Mg-values in the range —22.5 < Mb — 5 log /175 < M s (where M s = —15.12), in order to 
avoid small number statistics in each bin. 

We use the sets of distances r relative to different peculiar velocity field models, i.e., the 
multi-attractor model based on the whole Mark III data and on Mark III spiral data (Mark 
III*), Branchini & Plionis' (1996) cluster dipole model, the cluster dipole model as modified 
by Marinoni et al. (1998b) with the inclusion of a local model of the Virgocentric infall, the 
unperturbed Hubble flows in the CMB and Local Group (LG) frames. The number of galaxies 
[N ~ 5350) which fulfills the above-mentioned selection criteria changes according to the distance 
model adopted. 
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3 The galaxy luminosity function: results 

3. 1 The total luminosity function 

Table 1 reports the Schechter parameters obtained for sets of distance models. The errors 
given for a and M B are the projections onto their axes of the la confidence ellipse level, as 
derived from the \ 2 matrices of the least-square fit of Y(M#). 

The la formal errors in <p* include contributions both from galaxy clustering and from un- 
certainties in M B and a. We have checked that neglecting magnitude errors makes the slope a 
steeper by 0.03 and M* brighter by 0.1 mag, on average. 

The last columns of Table 1 contain the galaxy luminosity density pl = LQ(L)dL and 
the galaxy number density n together with the la errors. Letting the lower integration limit go 
to zero does not increase pL significantly. 

Fig. 1 show our results for the normalized LFs and the la error ellipses in the (a, Mb) plane, 
for galaxy distances predicted by the Mark III multi-attractor model and the cluster dipole 
model. Fig. 1 and the values of x 2 /dof given in Table 1 indicate that the Schechter function 
always gives a good fit to the LF of our sample. In the lowest luminosity bin there might be a 
hint for an excess of objects relative to the plotted fits for the Mark III multi-attractor models, 
but not for the cluster dipole model. 

Fig. 2 displays the observed galaxy number-distance histograms together with the curves 
which show the expected distribution for an uniform universe, for different distance models. 



TABLE 1 

The parameters of the Schechter LFs for different peculiar velocity field models. 



Model 


a 


M* - 5 log h 75 


4>* 


X 2 /dof 


PL 


n 








(10- 3 /if 5 Mpc- 3 ) 


(10 s 


' L Q h 75 Mpc 


" 3 ) (lO-^Mpc" 3 ) 


Multi-attractor (Mark III) 


-1.10 ±0.06 


-20.67 ± 0.08 


5.9 ±0.9 


0.52 


1.78 ± 0.28 


3.2 ±0.6 


Multi-attractor (Mark III*) 


-1.07 ±0.06 


-20.70 ± 0.08 


5.6 ±0.9 


0.42 


1.67 ±0.29 


2.9 ±0.5 


Cluster dipole model 


-1.15 ±0.05 


-20.63 ± 0.08 


5.9 ±0.9 


0.54 


1.75 ± 0.30 


3.7 ± 0.8 


Modified cluster dipole model 


-1.09 ±0.07 


-20.82 ± 0.08 


5.8 ±0.9 


0.68 


1.76 ± 0.33 


3.2 ±0.7 


Hubble Flow (CMB) 


-1.11 ±0.06 


-20.63 ± 0.07 


6.2 ±1.0 


0.49 


1.80 ±0.29 


3.5 ±0.7 


Hubble Flow (LG) 


-1.15 ±0.05 


-20.64 ± 0.08 


5.9 ±0.9 


0.24 


1.78 ± 0.31 


3.7 ± 0.8 
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M B -5 log h 75 

Fig. 1. — We plot the normalized differential LFs (together with la error bars) for galaxy distances derived 
from the multi-attractor model based on all Mark III data (left) and from the cluster dipole model data 
(right). 
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Fig. 2. — The plots show the observed galaxy number-distance histograms. They refer to redshift-distances 
obtained using the Hubble relation in the CMB frame (a) and in the LG frame (b), and to distances 
corrected according to the multi-attractor model fitted to the Mark III data (c), the multi-attractor 
model fitted to the Mark III* data (d), the cluster dipolc model (e), and the modified cluster dipolc 
model (f). The superimposed solid curves indicate the predictions given by the best Schechter function 
fits in the case of a homogeneous galaxy distribution. 
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The details of the large "peak and trough" fluctuations, which are seen both in wide-angle 
(e.g., the CfA2 sample) and narrow-angle surveys (e.g., the ESP sample by Zucca et al. 1997) 
are quite different in the various plots. Thus, the effects of peculiar motions have a large impact 
on the evaluation of the local galaxy density, especially on the smallest scales. On large scales, 
an overdensity at a distance of ~1500-2000 km s _1 (corresponding to the center of the Local 
Supercluster) and an underdensity around ~3500 km s" 1 are common to all histograms. 

All cases lead to similar values of (f>*, pl, n, and a. All cases give also similar M^-values, 
except for the modified cluster dipole model. Owing to its large Virgocentric infall, the modified 
cluster dipole model moves nearby galaxies quite close to us (see Fig. 2f). Moreover, it gives a 
M^-value typically brighter by ~0.1 - 0.2 mag than the norm, in agreement with the discussion 
by Efstathiou, Ellis & Peterson (1988) on the effect of a comparable Virgocentric infall on the 
LF of the RSA galaxy sample (Sandage & Tammann 1981). However, in recent years there is 
a growing evidence that the Virgo cluster is not a major source of the peculiar velocity field 
in the Local Supercluster. From the analysis of the magnitudes of bright cluster galaxies in a 
way which is free of assumptions about motions in the LS, Gudehus (1995) found no significant 
evidence of a Virgo infall of the LG. Therefore, in what follows we do not longer use the modified 
cluster dipole model, which is likely to be the least realistic model. 

To sum up, our calculations prove that differences in the peculiar velocity field models, though 
appreciably affecting the galaxy density field on galaxy-galaxy distance scale, do not have large 
effects on the optical LF of a shallow galaxy sample which covers a very large solid angle. The 
effects could be greater in samples restricted to fairly narrow solid angles (see, e.g., Davis & 
Huchra 1982, for the effect of a Virgocentric infall on the northern part of the CfAl sample). 

The values of pi and n tabulated in Table 1 refer to the population of cataloged galaxies 
of highish surface brightness. In general, magnitude-limited redshift surveys miss an increasing 
fraction of LSB objects at fainter magnitudes, which makes it difficult to determine the faint end 
of the LF. Among specific investigations designed to recover the LSB field population (e.g., the 
review by Impey & Bothun 1997), the large APM LSB redshift survey allowed the reconstruction 
of the LF of LSBs (Sprayberry et al. 1997), which shows that the LSB population (dominated 
by late- type spirals and irregulars), becomes appreciably numerous for magnitudes fainter than 
Mb ~-15 only. 

Hence, missing LSB objects do not affect appreciably the bright end of our LF. In the faint 
end they contribute strongly to the number density of galaxies (presumably also to cosmological 
mass density in baryons; e.g., Bristow Sz Phillipps 1994), but not too much to the luminosity 
density. Sprayberry et al. (1997) evaluated a contribution of 0.2 /175 L Q Mpc~ 3 (in the B 
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system) from the LF of the APM LSB galaxies, whereas Dalcanton et al. (1997a) derived 
PL = 0.1 hj^ Lq Mpc~ 3 (in the B system) from their small sample of late-type LSB galaxies 
identified in a large-area transit scan CCD data. 

With the addition of the ~10% light contribution given by Sprayberry et al. (1997), our 
mean value of pl = 2.0 ± O.3L /175 Mpc~ 3 implies a critical mass-to-light ratio (77^ = 
1.563-10 11 h 75 = 7>g ± 12 . 1Q 2 (M/L) Q to close the universe. The cosmological density 
parameter n = (M/L)/(M/L) C is ~0.2-0.3 if we set M/L ~ 200 h 75 ■ (M/L) , a typical virial 
mass-to-light ratio estimated from galaxy clusters in the blue (e.g., Bahcall, Lubin & Dorman 
1995; David, Jones & Forman 1995; Girardi et al. 1999). However, on larger scales, our analysis 
of cosmic flows (see Paper I) leads to a larger value of Qq, f^o ~0.5 (see also the reviews by 
Dekel, Burstein & White 1997 and Hamilton 1998). 



3.2 The type-specific luminosity functions 



Lumping together giants and dwarfs, we divide the NOG galaxies into five type bins, the 
ellipticals (E), lenticulars (SO), the Sa-Sb spirals (which cover the range between SO/a and Sb), 
the Sc-Sd spirals (which cover the range between Sbc and Sd), and the Sm-Im objects. We also 
consider the division into two type bins, the E-SO objects and the spirals and irregulars (S-Im). 
The morphological types are available for almost all the NOG galaxies. 

In Table 2 we give the results for the type-specific LFs (over the same ranges of distances and 
absolute magnitudes used for the whole sample), for different sets of distance models. In Fig. 3 
we plot the type-specific LFs and the total LF, in the specific case of Mark III distances. In Fig. 
4 We plot the lcr confidence ellipses for the parameters a and relative to the LFs shown in 
Fig. 3 for different morphological types. 

The morphological type-dependence of the LF is little sensitive to the distance models 
adopted. Assembling together the E and SO galaxies, we obtain a E-SO LF which remains 
almost flat, although the LF of the giant objects alone can be described by a Gaussian function 
(e.g., Binggeli et al. 1998; Biviano et al. 1995). More precisely, the LF of E galaxies appreciably 
declines towards faint luminosities (with a ~-0.5), whilst the LF of SO objects does not. 

The number of Sa-Sb types declines towards fainter magnitudes according to a slope (a ~ 
-0.6 - -0.8) which is considerably flatter than the norm for spirals, whilst the Sm-Im objects, 
which are on average, fainter than the norm, always exhibit a very steep slope (a ~-2.3 - -2.4). 
The slope a tends to get progressively steeper as one goes from early spirals to late spirals and 
irregulars. The value of MJ> is fainter for SO galaxies than for the other types. 
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TABLE 2 

The parameters of the morphological type-specific LFs. 



Model 


Sample 


Ngal 


a 


Mg — 5 log h 75 


0*(lO" 3 /i? 5 Mp C - 3 ) 


X 2 /dof 




E 


344 


-0.47 ±0.22 


-20.75 ±0.26 


0.46 ±0.12 


0.56 




SO 


596 


-1.17 ±0.20 


-20.44 ±0.26 


0.81 ±0.20 


0.37 




E-SO 


940 


-0.97 ±0.14 


-20.69 ±0.18 


1.03 ± 0.24 


0.51 


TV !T 1 i ■ ■ i i /"\ «" 1 TTT\ 

Multi-attractor (Mark III) 


Sa-Sb 


1521 


—0.62 ± 0.11 


—20.51 ± 0.12 


2.20 ± 0.46 


0.34 




Sc-Sd 


2240 


-0.89 ±0.10 


-20.39 ±0.11 


3.12 ± 0.59 


0.44 




Sm-Im 


619 


-2.41 ±0.28 


-21.11 ±0.72 


0.07 ±0.07 


0.68 




S-Im 


4380 


-1.10 ±0.07 


-20.63 ±0.09 


4.58 ± 0.73 


0.63 




E 


345 


-0.56 ± 0.22 


-20.71 ±0.26 


0.45 ±0.12 


0.39 




SO 


605 


-1.03 ±0.21 


-20.20 ±0.24 


1.03 ±0.26 


0.31 




E-SO 


950 


-1.03 ±0.14 


-20.67 ±0.19 


1.06 ±0.25 


0.64 


Cluster Dipole Model 


Sa-Sb 


1563 


-0.73 ± 0.11 


—20.48 ± 0.12 


2.24 ± 0.47 


0.49 




Sc-Sd 


2289 


-0.97 ± 0.09 


-20.35 ±0.12 


3.17 ± 0.60 


1.09 




Sm-Im 


597 


-2.45 ± 0.32 


-21.12 ±0.73 


0.07 ± 0.07 


0.64 




S-Im 


4449 


-1.17 ±0.08 


-20.60 ±0.09 


4.52 ± 0.72 


0.89 




E 


346 


-0.55 ±0.22 


-20.73 ±0.27 


0.45 ±0.11 


0.77 




SO 


589 


-1.08 ±0.21 


-20.18 ±0.24 


0.79 ± 0.20 


0.27 




E-SO 


935 


-1.12 ±0.13 


-20.72 ±0.19 


0.95 ±0.23 


0.74 


Hubble Flow in LG frame 


Sa-Sb 


1539 


-0.78 ±0.11 


-20.57 ±0.13 


2.10 ± 0.44 


0.57 




Sc-Sd 


2249 


-0.93 ± 0.09 


-20.33 ±0.11 


3.21 ±0.61 


0.37 




Sm-Im 


620 


-2.27 ±0.32 


-20.75 ±0.69 


0.15 ±0.15 


0.42 




S-Im 


4408 


-1.17 ±0.07 


-20.60 ±0.09 


4.62 ± 0.70 


0.45 



Except for the earliest (E) and latest (Sm-Im) types, there are no large differences between 
the shapes of the type-specific LFs. In particular, the shape of the E-SO LF does not differ 
significantly from that of the S-Im objects. 
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Fig. 3. — We plot the normalized LFs for all galaxy types and for different morphological types, adopting 
galaxy distances as given by the multi-attractor model fitted to the Mark III data. 
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Fig. 4. — We plot contours representing the 1 a deviation from the best-fitting values of a and Mg relative 
to the LFs shown in Fig. 3 for different morphological types. 
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3.3 The luminosity function of group galaxies 



With the purpose of investigating on the effect of the presence of galaxy systems on the LF 
determination, we recalculate the LFs after removing the members of the systems containing 
more than 5, 10, and 20 galaxies (e.g., removing 114, 42, and 10 systems for a total of 1437, 830, 
and 355 members, respectively, in the specific case of the Mark III multi-attractor model). We 
never find any significant difference in a and with respect to the values tabulated in Table 
1 for the whole sample. 

Furthermore, we calculate the LFs for galaxy members of groups and for field (non- grouped) 
galaxies, separately, according to the various distance models. Table 3 lists the results. We 
find no significant differences between the LF derived from the two subsamples, whenever, for 
constructing the sample of grouped galaxies, we take all systems or a subsample of systems 
dominated by groups with few members (e.g., all the systems with n > 5 members). For all sets 
of distance models, we detect some marginal systematic effects (at the ~ la confidence level) if 
we restrict ourselves to the richest systems containing n > 10, n > 15, n > 20 members. 

TABLE 3 



The Schechter LFs for 


FIELD 


AND GROUP 


GALAXIES. 




Model 


Sample 


Ngal 


a 


M* — 5 log h 75 


X 2 /dof 




field 


2652 


-1.19 ±0.10 


-20.59 ±0.12 


0.23 




groups 


2668 


-1.02 ±0.07 


-20.77 ±0.10 


0.58 


Multi-attractor (Mark III) 


groups(n > 10) 


963 


-1.21 ±0.11 


-20.99 ±0.18 


0.58 




groups(n > 20) 


394 


-1.28 ±0.18 


-21.00 ±0.31 


0.28 




field 


2735 


-1.24 ±0.09 


-20.57 ±0.12 


0.59 




groups 


2664 


-1.07 ±0.08 


-20.69 ±0.11 


0.38 


Cluster dipole model 


groups(n > 10) 


937 


-1.20 ±0.12 


-20.95 ±0.23 


0.62 




groups(n > 20) 


416 


-1.07 ±0.17 


-20.83 ±0.28 


0.49 




field 


2637 


-1.18 ±0.09 


-20.63 ±0.12 


0.25 




groups 


2706 


-1.15 ±0.07 


-20.70 ±0.11 


0.43 


Hubble Flow in LG frame 


groups(n > 10) 


952 


-1.35 ±0.11 


-20.97 ±0.19 


1.06 




groups(n > 20) 


395 


-1.37 ±0.17 


-21.05 ±0.44 


0.49 
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TABLE 4 

The parameters of the LFs from various samples. 



Sample 



Reference 



Ngal m u 



M min - 5 log h 75 



M* - 51ogfc 7 5 4>*{hhMpc~ 3 ) 



LCRS 
Autofib 
SSRS2 
ESP 
CS 

Durham/UKST Ratcliffe et al., 1998 



9100 mz = 15.5 



-18700 Gunn r~ 17.4 



Stromlo-APM Loveday et al., 1992 
CfA2 Marzke et al., 1994 

Lin et al., 1996 
Ellis et al., 1996 

Marzke & da Costa, 1997 ~3300 B(0) = 15.5 
Zucca et al., 1997 



Geller et al., 1997 



1700 R = 16.13 



-15.9 
-17.9 
-17.4 
-14.9 
-14.9 



-15.6 



-15.9 



-0.97 ±0.15 -20.42 ±0.13 (0.59 ± 0.07) ■ 10" 2 

-1.0 ±0.2 -20.2 ±0.3 (1.7 ±0.4) ■ 10~ 2 

-0.70 ±0.05 -20.11 ±0.02 (0.8 ± 0.04) • 10~ 2 

-1.16 ±0.05 -2O.22to.i2 ( 1 -04t ( o.i3) ' 10_2 

-1 - 16 to'.06 - 20.37 ±0.08 (0.46 ± 0.13) ■ 10" 2 



-13.3 -1.221°;* -20.53t°'og (0.8 ± 0.02) • 10" 



-1.17 ±0.19 -20.351°'^ (1.1 ±0.3) ■ 10" 2 
-1.04 ±0.08 -20.60 ±0.10 (0.72 ± 0.01) • 10" 2 



The LF of their members tend to show somewhat brighter values of Mg with respect to the 
field value, with no systematic shift of a. Although the average number of members of a group 
tends to decrease with increasing distance, the tendency is not so strong to affect significantly 
the above-mentioned results. 

4 Discussion 

4-1 Comparison with previous total galaxy luminosity functions 

In Table 4 we show the Schechter function parameters of some recent redshift surveys of 
field galaxies. All published values have been reduced to Hq = 75 km s" 1 Mpc" 1 . We quote 
the original photometric passband of the survey, but we transform the characteristic absolute 
magnitude M* to M|j through simple magnitude offsets; namely bj — r = +1.1, which is the 
mean rest-frame color of the Las Campanas redshif survey (LCRS) galaxies (Lin et al. 1996), 
bj — R = +1.3, which is the median rest-frame color of the Century Survey (CS) galaxies (e.g., 
Buta & Williams 1995; Geller et al. 1997), bj - m z = -0.45 (the shift of the blue Zwicky 
magnitudes mz comes from the comparison of galaxy number-magnitude counts in Shanks et 
al. 1984); B - bj = -0.3 as given by Efstathiou et al. 1988; B - B(0) = -0.3 (the zero point of 
the scale of the B(0) magnitudes agrees with the bj one within 0.2 mag, according to da Costa 
et al. 1994). The formal errors given in Table 4 do not include uncertainties related to color 
transformations. In the following, we quote the results given in the literature as scaled into the 
color and Hq we are using. 
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Fig. 5. — This figure shows a comparison between the NOG LF we calculate in the case of the multi- 
attractor model based on Mark III data, and previous LFs obtained in the literature from recent redshift 
surveys of low depth (a) and intermediate depth (b). 



Fig. 5 illustrates the comparison of our normalized LF (relative to the case of Mark IE-based 
distances) with those obtained from some low-depth (Fig. 5a) and intermediate-depth (Fig. 5b) 
redshift surveys of field galaxies. 

An accurate comparison between the LFs is impeded by the fact that applying simple constant 
magnitude offsets does neglect more complicated effects, which are due for example to galaxy 
colors, morphologies, and surface brightnesses. However, among the LFs there are certainly 
meaningful differences which are due to problems mentioned in §2. In any case, in the optical and 
near-infrared surveys mentioned in this section there is no indication of a systematic flattening of 
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the slope a, a simultaneous brightening of M B , and a progressive decrease of the normalization 
factor with brighter surface brightness limits of the survey, as expected in the case in which 
observational selection effects depending on brightness surface limits were the dominant effects 
(e.g., Dalcanton, Spergel & Summers 1997b). 

The M^-values given in Table 4 are always appreciably fainter (by ~0. 1-0.6 mag) than our 
values, which give a slower decline of the number of luminous galaxies with increasing luminosity. 

Our a-values well agree with the average results obtained from the shallow surveys, such as 
the Stromlo-APM, CfA2, SSRS2, CS, Durham/UKST surveys, and from intermediate-depth 
redshift surveys, such as the ESP and Autofib survey (restricted to 0.02 < z < 0.15). The only 
very discrepant value of a is that coming from LCRS. Its declining faint-end slope (a ~-0.7) 
can not be due to the fact that galaxies are selected in the red band, because the R-selected CS 
sample yields a much steeper slope. It can be partly related to the fact that a Schechter form 
does not fit very well the LCRS data around Mb ~-19.5 and in the faint end (see Bromley et 
al. 1998) and mostly to the fact that LCRS may be biased against observing low-luminosity 
galaxies, because of the survey selection criteria (Geller et al. 1997). 

In the intermediate luminosity interval (Mb ~-20) , our results point to a fairly low amplitude 
for the LF (see Fig. 5), in substantial agreement with the results coming from several shallow 
redshift surveys, especially with the Durham/UKST survey, which has yielded the LF which 
most closely resembles our results. 

It is interesting, albeit less easy, to compare our results with those coming fron the ORS. As a 
matter of fact, as warned by Santiago et al. (1996), it is difficult to properly compare the ORS 
LFs with other published LFs in the literature (including our results), because the ORS LFs are 
not corrected for Malmquist bias. Fig. 6 shows the comparison of our normalized LF (relative to 
the case of Hubble-flow distances in the LG frame) with the ORS LFs directly obtained from the 
UGC and ESO subsamples by Santiago et al. (1996), who parameterized the ORS LFs through 
a generalization of the Schechter form. For this comparison we apply no magnitude offset (since 
the ESO B system does not deviate too much from the standard RC3 B system; see e.g., Paturel, 
Bottinelli & Gouguenheim 1994), and no correction for differences in the mean galaxy density 
between the UGC and ESO subsamples. The NOG LF appears to be in satisfactory agreement 
with the ORS LFs at intermediate and faint magnitudes, whereas the bright end of the former 
seems to be shifted brightwards (like in the case of Fig. 6), suggesting magnitude offsets of 
~0.5 and ~0.7 mag, for the ESO and UGC subsamples respectively. Notably, the latter offset 
(which refers to Zwicky magnitudes) is much smaller than the amount (~ 1.2- 1.3 mag ) required 
to match the NOG LF and the bright end of the CfA2 LF (which also uses Zwicky photometry). 
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Fig. 6. — This figure shows a comparison between the NOG LF that we calculate in the case of redshift- 
distances obtained using the Hubble relation in the Local Group frame and the ORS LFs obtained directly 
from the ESO and UGC subsamples of the ORS. 

The anomalously high value of the CfA2 normalization is likely to be due to strong local 
inhomogeneities (e.g., the "Great Wall"). The quoted CS normalization, being the average 
value for the velocity range between 18 000 and 28 000 km s _1 , is also high, because in this 
range this survey slices through a portion of the Corona Borealis supercluster. Remarkably, both 
CfA2 and CS evidence that 4>* increases by a factor ~2 from nearby regions (cz < 5000 km s _1 ) 
to distant regions (cz ~ 10000 km s _1 and ~ 25000 km s _1 ), thus admitting a normalization 
consistent with our value in the nearby regions. 

The deeper redshift surveys (ESP, LCRS, Autofib) disagree in the amplitude of the LF at 
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intermediate and faint magnitudes. Our results are intermediate with respect to the low LCRS 
value and the high values of the ESP and Autofib samples. Zucca et al. (1997) stressed the 
possibility of reconciling the high ESP normalization with the low values given by other surveys, 
by noting that the nearby region (cz < 14000 km s" 1 ) encompassed by the ESP is underdense 
and that a recalculation of the LF for a subsample of bright ESP objects would gave results 
consistent with those derived from the Stromlo-APM survey. This view seems to be qualitatively 
supported by the preliminary results of the intermediate-depth 2dF redshift survey (whose 
galaxies are selected from the same catalog used in the last survey); also this LF hints at a 
greater normalization than the Stromlo-APM one (Maddox 1998). 

Noticeably, also the small-field deep near- infrared by Glazebrook et al. (1995b) and Cowie et 
al. (1996) disagree in the K-band LF normalization. The normalization found from the latter, 
larger survey is close to that coming from shallow near-infrared surveys (Mobasher, Sharpies & 
Ellis 1993; Gardner et al. 1997; Szokoly et al. 1998) and is few times smaller than that of the 
former survey. 

From the shallow IRAS 1.2 Jy survey conflicting results about the presence of a local void 
have emerged (cf Koranyi & Strauss 1997 and Springel & White 1998), whereas the deeper 
QDOT-IRAS survey suggests a moderate local void extending out to ~9000 km s" 1 (Saunders 
et al. 1990, Fig. 11). More reliable results may be expected from the larger PSCz survey (e.g., 
Canavezes et al. 1998) which is under study. 

To sum up, some optical redshift surveys (CfA2, CS, ESP) favor a high average normalization, 
but at the same time hint at a galaxy underdensity in the nearby regions. If we rely on the results 
coming from these surveys, we should conclude that the local universe is really underdense. But 
optical and IRAS redshift surveys disagree on the extent, depth, and even real existence of a 
local underdensity (see §4.4 for further discussions on this point). 

4-2 Comparison with previous type-specific galaxy luminosity functions 

Our morphological-type dependence of the slope of the optical galaxy LF is in substantial 
agreement with the results obtained by Marzke et al. (1994b, 1998) from the CfA2 and the 
SSRS2 samples. The number of Sm-Im objects rises towards the faint end even steeper in our 
sample than in those two samples (where a ~-1.9). Moreover, the difference between the a- values 
of E and SO galaxies (and correspondingly the shortage of low-luminosity E systems) is more 
marked in our sample than in the CfA2 sample. The difference we find is fully consistent with 
claims that in the field the transition region between giants and dwarfs is scarcely populated in 
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the case of ellipticals and better populated in the case of lenticulars, because galaxy bulges tend 
to weaken with decreasing luminosities and become dominated by disk or lens-like components 
(see, e.g., the RSA survey of dwarf galaxies by Vader & Chaboyer 1994). 

Nevertheless, the dependence of the slope on the morphological type across the whole Hubble 
morphological sequence is not so striking in our sample as in the Stromlo-APM sample, where the 
very steep declining of early-type galaxies at the faint end (Loveday et al. 1992) can be ascribed 
to the incompleteness of early types at high redshift and/or faint magnitudes (as discussed by 
Marzke et al. 1994b, and Zucca, Pozzetti & Zamorani 1994). 

Our analysis establishes that the M^-value of E galaxies is not unusually bright, supporting 
Marzke et al.'s (1994b) guess that an uncorrected dependence of the Zwicky magnitude scale on 
surface brightness (e.g., Ichikawa & Fukugita 1992) could have bias the magnitudes of the CfA2 
ellipticals. 

At middle and faint magnitudes (Mb >-20), our spiral LF has an amplitude closer to those 
of the CfA2 and SSRS2 spiral LFs than to the lower one of the Stromlo-APM late-type LF. At 
the same magnitudes, the amplitude of our E-SO LF is close to that of the SSRS2 E-SO LF. 

Interestingly, the spectral type-dependence of the galaxy luminosity function as derived from 
the LCRS ( Bromley et al. 1998 a, b) and the Autofib redshift survey (Heyl et al. 1997) shows a 
more progressive steepening of the slope from early types to late types than the morphological- 
type dependence derived in this paper. This is due to the fact that galaxy typing based on the 
visible morphological appearance is not very tightly correlated with galaxy classification based 
on spectral schemes (e.g., Kennicutt 1992; Connolly et al. 1995; Zaritsky, Zabludoff & Willick 
1995; Katgert et al. 1999) because the latter method is very sensitive to the degree of star 
formation activity. 

Our study lends further support to the fact that blue, gas-rich, late-type spirals and irregulars 
constitute a steeply rising component of the field galaxy LF. Since these galaxy types dominate 
in the field at very low luminosities (e.g., Ferguson & Binggeli 1994; Sprayberry et al. 1997), our 
finding leads us to envisage an upturn of the field LF in the faint end (though not brightwards 
Mb ~-15). This is consistent with the growing evidence that the field galaxy LF shows an 
upturn from the Schechter form at faint magnitudes (where magnitude-limited samples become 
severely censored by surface brightness selection effects). There is straightforward evidence 
for this faint-end steepening in the CfA2 and ESP data (through an excess of low-luminosity 
galaxies above the extrapolated Schechter function) and in the counts of the statistical excess of 
faint APM galaxies seen in projection around Stromlo-APM redshift survey galaxies (Loveday 
1997). 
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4-3 Comparison with previous luminosity functions of group and cluster galaxies 



The most recent and extensive analyses of the cluster galaxy LF (Lumsden et al. 1997; Valotto 
et al. 1997; Gaidos 1997; Rauzy, Adami & Mazure 1998) concur to M^j-values brighter by a 
few tenths of magnitudes than those of the field galaxy LFs mentioned in §4.1 (including our 
LF and the LFs relying on the same photometric system), whereas there is no agreement on 
whether the LF slope a of members of (poor and/or rich) clusters significantly differs than that 
of field galaxy LF (see the conflicting results reached by the above-mentioned authors and also 
by Lopez-Cruz et al. 1997). 

Notably, that tendency is consistent with the environmental dependence of the LF we recognize 
in §3.3. We find that this dependence is so much weak (at least in the systems of moderate 
richness which are typical for our sample) that the presence of galaxy systems does not affect 
significantly the field galaxy LF. This is at variance with some pronounced cluster effects (i.e., a 
steepening of the slope a in samples containing clusters) reported by Marzke & da Costa (1997) 
for a small subsample drawn from the SSRS2. 

Our results are in line with the claims of some luminosity segregation with density for several 
galaxy samples (e.g., Park et al. 1994, for the CfA2 galaxies; Giuricin et al. 1995, for the 
Local Supercluster galaxies; Loveday et al. 1995 for the Stromlo-APM galaxies, Willmer, da 
Costa &: Pellegrini 1998 for the SSRS2 galaxies). This segregation was found to be weaker 
than morphological segregation with density by Iovino et al. (1993) for the Perseus-Pisces 
supercluster galaxies and by Dominguez-Tenreiro, Gomez-Flechoso & Martinez (1994) for the 
CfAl galaxies, whereas it was considered significant only for early-type luminous galaxies by 
Hasegawa & Umemura (1993). 

Among recent works in which environmental effects on the LF are searched for through a 
specific comparison between group and field galaxy LFs, an approach similar to ours was followed 
by Ramella et al. (1999) and Bromley et al. (1988a,b), who divided their respective ESP and 
LCRS population into field galaxies and galaxy members of groups (or rich groups), identified 
by means of the friends of friends algorithm. 

The former authors found that, compared to the LF of field galaxies, the LF of group galaxies 
exhibits a brighter M* (by a few tenths of magnitude), which brightens with the increasing 
richness of groups, whereas the slope a does not change significantly from the field to the 
groups. Ramella et al.'s (1999) results on changes of M* are quantitatively compatible with our 
findings, although the authors claimed a larger (2 a) statistical significance. 

On the other hand, the latter authors detected a steepening of the LF slope for galaxies 
located in regions of higher density. This is caused by a strong environmental effect (in the same 
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sense) on the LF of early-type galaxies (there is no effect on late types). We must be wary that 
their steepness-density relation (which we do not find) could reflect some property of the survey 
(which shows a north-south discrepancy in total and type-dependent LFs) rather than a true 
environmental effect. 

4-4 Relation of the galaxy luminosity function to galaxy counts 

The following arguments bearing on the observed galaxy number counts provide a better 
evidence for an underdensity of the nearby universe than those presented in §4.1. As discussed 
by several authors (e.g., Driver & Phillipps 1996, and the reviews by Shanks 1989, and Ellis 
1997), in the absence of appreciable galaxy evolution, a low normalization of the local galaxy 
LF, as derived by us, is inconsistent with the high counts which are generally observed already 
at intermediate optical and near-infrared magnitudes (i.e., roughly at 18 < B < 20 mag, 13 < 
K < 15 mag). The observed counts at these magnitudes would imply a LF normalization 
greater by a factor ~1.5 than that of our LF. A high normalization is customarily assumed in 
most models of faint galaxy counts, which are generally normalized to the observations in the 
magnitude intervals quoted above rather than locally. Workers (e.g., Koo, Gronwall & Bruzual 
1993; Gronwall & Koo 1995), who claimed that the high bright-end counts can be reasonably 
fitted by no-evolution or mild evolution models, adopted the view that the local LF is ill-defined. 

We deem it more reasonable to explain the underprediction of galaxy counts by conceiving a 
local underdensity rather than a rapid evolution (Maddox et al. 1990) of the bulk of the galaxy 
population at the corresponding low redshifts of z ~0.1 - 0.2 because on theoretical grounds 
little galaxy evolution is expected up to these redshifts and because very deep redshift surveys 
which have begun measuring the galaxy LF out to z ~0.5 - 1 (e.g., Ellis 1997) reveal little 
evolution in the LF of most of galaxies, except for the late-type low-luminosity population. 

Consistently, no-evolution models based on our local LF , but with a normalization roughly 1.5 
times our value, would satisfactorily predict the counts of E/SO and Sabc spirals, but would still 
underpredict those of late-type spirals and irregulars, for the morphologically classified sample 
of the faint field galaxies from the HST Medium Deep Survey (Driver, Windhorst & Griffiths 
1995; Glazebrook et al. 1995a; Abraham et al. 1996) and the deep HST fields (Odewahn et al. 
1996). 

As for the compatibility of the observed bright end counts with no-evolution or passive evo- 
lution models, recent observations disclaimed the very steep slope of the APM counts reported 
in the 16 < bj < 19 mag range by Maddox et al. (1990) — who explained their counts by 
invoking a large amount of evolution at low redshift — and attributed the steepness of the APM 
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counts to systematic errors in the non-linear APM photometry. The same effect may account for 
the low counts of Heydon-Dumbleton, Collins & MacGillivray (1989) and Zucca et al. (1997) at 
bright COSMOS-frj magnitudes. The flatter slopes which fit the bright-end blue counts of recent 
studies (Weir, Djorgovski & Fayyad 1995; Gardner et al. 1996, Huang et al. 1997, Bertin k, 
Dennefeld 1997) do not largely depart from the expectations of no-evolution models, but always 
require a higher normalization than our local one. 

A more convincing evidence for such a departure comes from the bright-end K-band number 
counts. At bright magnitudes in the K-band (10 < K < 17 mag), where the K-correction is 
known better than in the optical band and is only a weak function of morphological type, the 
number counts delineate a slope which is too steep to be consistent with simple no-evolution 
models (Huang et al. 1997). Detailed models of these counts have explained this effect through 

a large region of local underdensity of galaxies (by a factor 1.7 2.4), which extends out to 

z ~0.1 - 0.2, in the case of no-evolution, or through a combination of some local underdensity 
and a mild evolution (Phillips & Turner 1998). 

The initial steep rise of the radio source counts between fluxes of ~10 Jy and ~1 Jy (in the 
1.4 -8.4 GHz frequency range) (e.g., Windhorst et al. 1993) also may suggest the presence of a 
local underdensity of nearby radiogalaxies, since their strong cosmological evolution is thought 
not to start until z ~0.3 (e.g., Condon 1989). 

The existence of a large-scale local underdensity implies that local measurements underesti- 
mate (overestimate) the true value of fto (Hq). The local value of Hq is expected to be higher 
than the global value (Turner, Cen & Ostriker 1992) by as much as ~20% for f^o =1 an d ~10% 
for rio = 0.2, in the case of an underdensity by a factor of 1.5 (see also Shi, Widrow & Dursi 
1996) . A moderate density contrast (of this order of magnitude or smaller) is lower than the value 
favored by the above-mentioned models of bright-end K-band galaxy number counts (Phillips 
k, Turner 1998), but it is more consistent with the following observational constraints relative 
to Hubble flow perturbations. As a matter of fact, from the peculiar velocities of supernovae 
Zehavi et al. (1998) found a ~7% deviation from the Hubble law consistent with a void of ~20% 
underdensity surrounded by a dense wall at 7000 km s _1 (which roughly coincides with the local 
Great Walls). This is compatible with the corresponding upper limit (~7%) obtained from the 
brightest cluster galaxy Hubble diagram (Lauer k Postman 1992; Lauer et al. 1998). Differences 
between the local and global values of Hq as large as ~10% on scales of ~30,000 km s _1 are 
still fully compatible with the limits sets by the measurements of cosmological distances through 
high redshift supernovae (e.g., Kim et al. 1997), but become only marginally compatible with 
the limits imposed by the current knowledge of cosmic microwave background anisotropics (e.g., 
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Wang, Spergel & Turner 1998). 



5 The selection function 



The NOG sample, being magnitude-limited, becomes sparser at larger distances due to the 
increasing loss of galaxies caused by the apparent magnitude cutoff. This effect is quantified by 
the selection function S(r) which expresses the fraction of galaxies that are expected to satisfy 
the sample's selection criterion. The incompletion factor F(r), which expresses the number of 
galaxies that should have been cataloged for each object present in the sample at a given distance 
r, is related to S(r) by the expression F(r) = 1/S(r). 

The selection function is given by 



where 3> is the convolved Schechter LF, / = 0.8 is the average completeness level, L m i n (r) is the 
minimum luminosity necessary for a galaxy at distance r (in Mpc) to make it into the sample, and 
the integral is cut off at the lower limit of L s = L m i n (r s ) corresponding to M s = —15.12+5 log /175 
(see §2.2). Through Monte Carlo simulations Santiago et al. (1996) 

demonstrated that an unbiased estimate of the galaxy density field can be obtained in the 
presence of random errors on magnitudes if one uses the convolved LF in eq. [LT|. 

Fig. 7 shows the incompletion factor F as a function of the distance r (in km s" 1 ) for different 
velocity field models. The incompletion factor F is always a smaller in the case of the Mark III 
multi-attractor model than in the case of the cluster dipole model, but F appears to be little 
sensitive to the adopted models. 

The dependence of the selection function on redshift-space distortions has been addressed 
analytically by Hamilton (1998) in the case of linear theory. His work is relevant in order to 
understand the results given in Fig. 7, and then it is mentioned here. 

Denoting the redshift-space coordinate by s, the real-space (F r (r)) and redshift-space (F s (s)) 
incompletion functions in the LG frame are related as follows (see eq. (4.85) in Hamilton 1998): 



where r max is the depth of the sample. Note that averages are not ensemble averages, in which 
case the real- and redshift-space F would be equal, but are performed over the actual galaxy 
sample. 




(11) 




(12) 
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Fig. 7. — We plot the incompletion factor F as a function of galaxy distance r (expressed in km/s) for 
different peculiar velocity field models. 



The second term inside the parenthesis of the right hand side of eq. 12 averages to zero for 
a full-sky catalog; this highlights a further advantage of the full-sky coverage. The first term 
in the same parenthesis is the average, at the depth r, of the line-of-sight peculiar velocity v 
over r. Then, the largest contribution to the differences between the real- and redshift-space 
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incompletion functions comes from the very nearby universe, both because of the 1/r term and 
because the velocity v is less affected by local fluctuations if the average volume is larger. 



With eq |12| , considering that the logarithmic derivative is positive for common flux-limited 
samples, it is possible to understand the behaviors of the different models. As a matter of fact, 
any model which gives a local overdensity has negative (v/r) r (converging velocities). This is 
the case for the Mark III multi-attractor model, which has a nearby attractor (the Virgo cluster) 
and two other attractors (the Great Attractor and the Perseus-Pisces supercluster) at distances 
of ~4000 - 5000 km/sec (the distant Shapley concentration gives a significant contribution to 
the LG motion, but covers a small solid angle). As a consequence, galaxies are brighter (more 
distant) on average in real space, and the corresponding incompletion factor is consequently 
smaller. 

On the other hand, the cluster dipole model gives a much smoother version of the local 
velocity field. As a consequence, the local overdensity is low and the incompletion function 
remains similar to the LG one. 

6 Summary and conclusions 

We evaluate the field galaxy LF and the related selection function of the nearly complete, 
magnitude-limited (B< 14 mag) and volume-limited (cz <5500 km s _1 ) NOG (Nearby Optical 
Galaxy) sample. Our derivation of the LF is performed using the locations of field and grouped 
galaxies in real distance space, since we adopt sets of galaxy distances based on different peculiar 
velocity field models. Our main results can be summarized as follows: 

i) The local field galaxy LF is well described by a Schechter function with a ~-l.l, <£* ~ 
0.006 Mpc~ 3 , M* B ~-20.6 down to M B = -15.1 (H = 75 km s' 1 Mpc' 1 ). The exact values 
of the Schechter parameters little depend on the adopted velocity field models, since peculiar 
motion effects are of the order of statistical errors (at most, they cause variations of 0.08 (~ la 
error) in a and 0.2 mag (~ 2a error) in Mg). Thus, the local field galaxy LF, as well as the 
intimately related selection function of the NOG sample, are little sensitive to peculiar motions 
effects. On the other hand, the details of peculiar velocity models have a quite large impact on 
the local galaxy density on the smallest scales. 

ii) Compared to previous local field galaxy LFs, our LF is characterized by a normal (mod- 
erately steep) slope a and a M^-value which is, on average, brighter by 0.4 mag than the 
norm. The last result is explained from the fact that the photometric measures of the NOG 
galaxies, being total B magnitudes corrected for Galactic absorption, K-dimming, and internal 
absorption, better represent the true galaxy light. 
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iii) Our morphological-type dependence of the galaxy LF is in much closer agreement with 
the results obtained from the CfA2 and SSRS2 samples (Marzke et al. 1994b, 1998) than with 
those derived from the Stromlo-APM sample (Loveday et al. 1992). The dependence of the 
galaxy LF from the morphological type appreciably differs from its dependence on the spectral 
classification as given in the literature. 

The LF shape of E-SO galaxies does not differ significantly from the that of spirals and irreg- 
ulars. On the other hand, the E galaxies clearly decrease in number towards low luminosities, 
whereas the number of late-type spirals and irregulars show a very steep rise towards the faint 
end (with a ~-2.3 - -2.4). This behaviour suggests a similarly fast increase in the number of 
low-luminosity galaxies in the faint end unexplored in the NOG sample (at Mb > —15). 

iv) The presence of galaxy systems in the NOG sample does not affect significantly the shape 
of the the field galaxy LF. As a matter of fact, environmental effects on the total LF are proved 
to be marginal in the NOG sample. However, the total LF of the galaxy members of the richest 
systems tends to show a somewhat brighter value of M B (at the la level), which is in line with 
several earlier claims of galaxy luminosity segregation with density (in particular with Ramella 
et al.'s (1999) results). 

This effect is unrelated to an increased proportion of early-type galaxies in galaxy systems. A 
luminosity segregation with density could be related to the processes of biased galaxy formation 
(e.g., Kaiser 1984; Rees 1985) and/or could be induced by late evolutionary processes — such as 
tidal shocks (e.g., Nogushi & Ishibashi 1986; Sanders et al. 1988), ram pressure effects (Dressier 
& Gunn 1983), and merging effects (e.g., Barnes & Hernquist 1991) — which can make galaxies 
more prone to starburst activity whenever they are located in denser environments (e.g., Maia 
et al. 1994), at least in cases of not extreme environments (Hashimoto et al. 1998, Balogh et 
al. 1998). 

v) Our results supports the fact the normalization of the local galaxy LF is relatively low. 
Hence, our local galaxy LF underpredicts many observed galaxy number counts (including those 
of E/SO galaxies) at relatively bright optical and near- infrared magnitudes (where little galaxy 
evolution is in general expected and observed). If we do not want to raise serious doubts on 
the basic aspects of our understanding of galaxy evolution and spectral energy distribution, the 
most plausible explanation for this underprediction is that the nearby volume we investigate is 
underdense in galaxies (by a factor of ~1.5). However, we are aware of the systematic uncertain- 
ties which in general tend to bias low the estimate of the LF normalization factor (e.g., Willmer 
1997) and of the appreciable sensitivity of this factor to density anomalies within local volumes. 

Notably, uncertainties related to the choice of the global LF normalization have a much larger 



32 



impact on the census of light in the universe (hence, on the determination of the critical mass- 
to-light ratio) than uncertainties in the light contribution of uncataloged LSB galaxies (see 
Dalcanton, Spergel & Summers 1997b for a dissenting view). 

In conclusion, large observational efforts (e.g., redshift surveys) are still needed to reliably 
establish the existence, density contrast, and size of a possible local void, which, as suggested 
by some authors, can be larger than the ~5000 km s _1 radius region examined in this paper. 
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